Search for: All records

Abstract Volcanic evolution in ocean island settings is often controlled by variations in the chemistry and volumetric flux of magma from an underlying mantle plume. In locations such as Hawaiʻi or Réunion, this results in predictable variations in magma chemistry, the rate of volcanic activity, and the depth of magma storage with volcanic age and/or distance from the centre of plume upwelling. These systems, however, represent outliers in global plume volcanism due to their high buoyancy flux, frequent eruptions, and large distance from any plate boundary. Most mantle plumes display clear interaction with nearby plate boundaries, influencing the dynamics of solid plume material in the upper mantle and the distribution of melt across regions of active volcanism. Yet, the influence of plume–ridge interaction and plume–ridge distance on the structure, characteristics, and evolution of magma storage beneath ocean island volcanoes remains under constrained. In this study, we consider the evolution of magmatic systems in the Galápagos Archipelago, a region of mantle plume volcanism located 150–250 km south of the Galápagos Spreading Centre (GSC), focusing on the depth of magma storage during the eastward transport of volcanic systems away from the centre of plume upwelling. Geochemical analysis of gabbro xenoliths from Isla Floreana in the southeastern Galápagos suggest that they formed at ~2–2.5 Ma, when the island was located close to the centre of plume upwelling. These nodules, therefore, provide rare insights into the evolution of volcanic systems in the Galápagos Archipelago, tracking variations in the magma system architecture as the Nazca plate carried Isla Floreana eastwards, away from the plume centre. Mineral thermobarometry, thermodynamic modelling, and CO2 fluid inclusion barometry reveal that Isla Floreana’s plume-proximal stage of volcanic activity—recorded in the gabbro xenoliths—was characterized by the presence of high-pressure magma storage (>25 km), below the base of the crust. In fact, we find no petrological evidence that sustained, crustal-level magma storage ever occurred beneath Isla Floreana. Our results contrast with the characteristics of volcanic systems in the western Galápagos above the current centre of plume upwelling, where mid-crust magma storage has been identified. We propose that this change in magmatic architecture of plume-proximal volcanic centres in the Galápagos—from high-pressure mantle storage at 2.5 Ma to mid-crustal storage at the present day—is controlled by the variations in plume–ridge distance. Owing to the northward migration of the GSC, the distance separating the plume stem and GSC is not constant, and was likely <100 km at 2.5 Ma, significantly less than the current plume–ridge distance of 150–250 km. We propose that smaller plume–ridge distances result in greater diversion of plume-material to the GSC, ‘starving’ the eastern Galápagos islands of magma during their initial formation and restricting the ability for these systems to develop long-lived crustal magma reservoirs. 

Mineral and fluid inclusion data collected from a suite a gabbro xenoliths found on the island of Floreana in the southern Galapagos. Mineral chemistry is primarily determined by SEM-based EDS and WDS analysis with LA ICP-MS used to determine the clinopyroxene trace element contents. Microthermometry and Raman analysis are used to assess the composition and CO2 density of clinopyroxene and plagioclase hosted fluid inclusions. 

Abstract Rapid-response petrological monitoring is a major advance for volcano observatories, allowing them to build and validate models of plumbing systems that supply eruptions in near-real time. The depth of magma storage has recently been identified as high-priority information for volcanic observatories, yet this information is not currently obtainable via petrological monitoring methods on timescales relevant to eruption response. Fluid inclusion barometry (using micro-thermometry or Raman spectroscopy) is a well-established petrological method to estimate magma storage depths and has been proposed to have potential as a rapid-response monitoring tool, although this has not been formally demonstrated. To address this deficiency, we performed a near-real-time rapid-response simulation for the September 2023 eruption of Kīlauea, Hawaiʻi. We show that Raman-based fluid inclusion barometry can robustly determine reservoir depths within a day of receiving samples—a transformative timescale that has not previously been achieved by petrological methods. Fluid inclusion barometry using micro-thermometric techniques has typically been limited to systems with relatively deep magma storage (>0.4 g/cm3 i.e.  > 7 km) where measurements of CO2 density are easy and accurate because the CO2 fluid homogenizes into the liquid phase. Improvements of the accuracy of Raman spectroscopy measurements of fluids with low CO2 density over the past couple of decades has enabled measurements of fluid inclusions from shallower magmatic systems. However, one caveat of examining shallower systems is that the fraction of H2O in the fluid may be too high to reliably convert CO2 density to pressure. To test the global applicability of rapid response fluid inclusion barometry, we compiled a global melt inclusion dataset (>4000 samples) and calculate the fluid composition at the point of vapor saturation ($${\mathrm{X}}_{{\mathrm{H}}_2\mathrm{O}}$$). We show that fluid inclusions in crystal hosts from mafic compositions (<57 wt. % SiO2)—likely representative of magmas recharging many volcanic systems worldwide—trap fluids with $${\mathrm{X}}_{{\mathrm{H}}_2\mathrm{O}}$$ low enough to make fluid inclusion barometry useful at many of the world’s most active and hazardous mafic volcanic systems (e.g. Iceland, Hawaiʻi, Galápagos Islands, East African Rift, Réunion, Canary Islands, Azores, Cabo Verde). 

We present DiadFit—an open-source Python3 tool for efficient processing of Raman spectroscopy data collected from fluid inclusions, melt inclusions and silicate melts. DiadFit is optimized to fit the characteristic peaks from CO2 fluids (Fermi diads, hot bands, 13C), gas species such as SO2, N2, solid precipitates (e.g. carbonates), and Ne emission lines with easily tweakable background positions and peak shapes. DiadFit's peak fitting functions are used as part of a number of workflows optimized for quantification of CO2 in melt inclusion vapour bubbles and fluid inclusions. DiadFit can also convert between temperature, pressure and density using various CO2 and CO2-H2O equations of state (EOS), allowing calculation of fluid inclusion pressures (and depths in the crust), conversion of homogenization temperatures from microthermometry to CO2 density, and propagation of uncertainties associated with EOS calculations using Monte Carlo methods. There are also functions to quantify the area ratio of the silicate vs. H2O region of spectra collected on silicate glasses to determine H2O contents in glasses and melt inclusions. Documentation and worked examples are available (https://bit.ly/DiadFitRTD, https://bit.ly/DiadFitYouTube). 

Raman spectroscopy has become the tool of choice for analyzing fluid inclusions and melt inclusion (MI) vapor bubbles as it allows the density of CO2-rich fluids to be quantified. Measurements are often made at ambient temperature (Tamb ~18-25 °C), resulting in reported bulk densities between 0.2 and 0.7 g/mL despite that single-phase CO2 under these conditions is thermodynamically unstable and instead consists of a liquid (~0.7 g/mL), and a vapor phase (~0.2 g/mL). Here, we present results from experiments conducted at Tamb and 37 °C (above the CO2 critical temperature) on 14 natural CO2-rich MI bubbles from Mount Morning, Antarctica. Here, we show that at Tamb, laser power strongly affects the CO2 Raman spectrum of MI bubbles with bulk densities within the miscibility gap. High-power laser heating and low spectral resolution explain why published measurements have reported such bulk densities at Tamb even when using an instrument-specific calibration. 

Major and trace element abundances, including highly siderophile elements, and 187Os and 182W isotopic compositions were determined for ca. 89 Ma mafic and ultramafic rocks from the islands of Gorgona (Colombia) and Curaçao (Dutch Caribbean). The volcanic systems of both islands were likely associated with a mantle plume that generated the Caribbean Large Igneous Provence. The major and lithophile trace element characteristics of the rocks examined are consistent with the results of prior studies, and indicate derivation from both a chemically highly-depleted mantle component, and an enriched, or less highly-depleted mantle component. Highly siderophile element abundances for these rocks are generally similar to rocks with comparable MgO globally, indicating that the major source components were not substantially enriched or depleted in these elements. Rhenium-Os isotopic systematics of most rocks of both islands indicate derivation from a mantle source with an initial 187Os/188Os ratio between that of the contemporaneous average depleted mid-ocean ridge mantle and bulk silicate Earth. The composition may reflect either an average lower mantle signature, or global-scale Os isotopic heterogeneity in the upper mantle. Some of the basalts, as well as two of the komatiites, are characterized by calculated initial 187Os/188Os ratios 10-15% higher than the chondritic reference. These more radiogenic Os isotopic compositions do not correlate with major or trace element systematics, and indicate a mantle source component that was most likely produced by either sulfide metasomatism or ancient Re/Os fractionation. Tungsten-182 isotopic compositions measured for rocks from both islands are characterized by variable 182W values ranging from modern bulk silicate Earth-like to strongly negative values. The 182W values do not correlate with major/trace element abundances or initial 187Os/188Os compositions. As with some modern ocean island basalt systems, however, the lowest 182W value (-53) measured, for a Gorgona olivine gabbro, corresponds with the highest 3He/4He previously measured from the suite (15.8 R/RA). Given the lack of correlation with other chemical/isotopic compositions, the mantle component characterized by negative 182W and possibly high 3He/4He is most parsimoniously explained to have formed as a result of isotopic equilibration between the mantle and core at the core-mantle boundary. 

Constraining the volatile content of magmas is critical to our understanding of eruptive processes and their deep Earth cycling essential to planetary habitability [R. Dasgupta, M. M. Hirschmann, Earth Planet. Sci. Lett. 298 , 1 (2010)]. Yet, much of the work thus far on magmatic volatiles has been dedicated to understanding their cycling through subduction zones. Further, studies of intraplate mafic volcanism have disproportionately focused on Hawaii [P. E. Wieser et al., Geochem. Geophys. Geosyst. 22 , e2020GC009364 (2021)], making assessments of the overall role of intraplate volcanoes in the global volatile cycles a challenge. Additionally, while mafic volcanoes are the most common landform on Earth and the Solar System [C. A. Wood, J. Volcanol. Geotherm. Res. 7 , 387–413 (1980)], they tend to be overlooked in favor of silicic volcanoes when it comes to their potential for explosivity. Here, we report primitive (olivine-hosted, with host Magnesium number – Mg# 78 to 88%) melt inclusion (MI) data from Fogo volcano, Cabo Verde, that suggest that oceanic intraplate silica-undersaturated explosive eruptions sample volatile-rich sources. Primitive MI (melt Mg# 70 to 71%) data suggest that these melts are oxidized (NiNiO to NiNiO+1) and very high in volatiles (up to 2 wt% CO 2 , 2.8 wt% H 2 O, 6,000 ppm S, 1,900 ppm F, and 1,100 ppm Cl) making Fogo a global endmember. Storage depths calculated from these high volatile contents also imply that magma storage at Fogo occurs at mantle depths (~20 to 30 km) and that these eruptions are fed from the mantle. Our results suggest that oceanic intraplate mafic eruptions are sustained from the mantle by high volatile concentrations inherited from their source and that deep CO 2 exsolution (here up to ~800 MPa) drives their ascent and explosivity.